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GENERALIZATION OF BOUNDARY-LAYER MOMENTUM-INTEGRAL EQUATIONS TO 


THREE-DIMENSIONAL FLOWS INCLUDING THOSE OF ROTATING SYSTEM! 


By Artur MAGER 


SUMMARY 


The Navier-Stokes equations of motion and the equation of 
continuity are transformed so as to apply to an orthogonal 
curvilinear coordinate system rotating with a uniform angular 
velocity about an arbitrary axis in space. A usual simplifica- 
tion of these equations as consistent with the accepted boundary- 
layer theory and an integration of these equations through the 
boundary layer result in boundary-layer momentum-integral 
equations for three-dimensional flows that are applicable to 
either rotating or nonrotating fluid boundaries. 

These equations are simplified and an approximate solution 
in closed integral form is obtained for a generalized boundary- 
layer momentum-loss thickness and flow deflection at the wall 
in the turbulent case. ; 

sl numerical evaluation of this solution carried out for data 
obtained in a curving nonrotating duct shows a fair quantita- 
tire agreement with the measured values. 

The form in which the equations are presented is readily 
adaptable to cases of steady, three-dimensional, incompressible 
boundary-layer flow like that over curved ducts or yawed wings; 
and it also may be used to describe the boundary-layer flow over 
various rotating surfaces, thus applying to turbumachinery, 
propellers, and helicopter blades. 


INTRODUCTION 


The development of the boundary layer on the various 
parts of turbomachinery (compressors and turbines), heli- 
copter blades, propellers, and in curved ducts is influenced by 
centrifugal and Coriolis forces in addition to the pressure and 
viscous forces. As a result of these forces, the flow in the 
boundary layer not only has the characteristic velocity 


deficiency but also has, because of this velocity deficiency, 


direction different from that of the flow outside the boundary 
layer. Thus the behavior of the boundary layer in three- 
‘dimensional flow may be quite unlike the behavior in two- 
dimensional flow. The effect of these additional forces on the 
boundary layer has been realized for some time and the 
observed discrepancies in the boundary-layer behavior have 
usually been explained only in a qualitative manner as, for 
example, in references 1 to 4. 

The literature concerning the theoretical aspect of the 
three-dimensional boundary-layer flow is meager. For the 
laminar case most of the published work has been carried out 
in connection with the yawed wing (references 5 to 7). For 
the turbulent case, although a number of researchers have 
established the general form of the differential equations 


‘out along a streamline of the potential flow. 





applicable, no actual solutions of these equations have been 
obtained. Tetervin, for example, presents boundary-layer 
momentum-integral equations in three dimensions for a fluid 
of variable density and viscosity (reference 8). Gruschwitz 
establishes the momentum-integral equations for boundary- 
layer flow along an arbitrarily curved streamline in reference 9. 
Burgers gives the differential equations on the develop- 
ment of boundary layers in the case of axially symmetric 


flows having a rotational component (reference 10). Prandtl, 


in addition to presenting a form of three-dimensional 
momentum-integral equations, suggests the general procedure 
that could be followed to obtain a solution (reference 11). 
Experimental data are similarly lacking. In spite of con- 
siderable literature search, only the data of Gruschwitz 
(reference 9) for a curved duct and the data of Kuethe, 
McKee, and Curry (reference 12) for a yawed wing were 
found. : 

As a result of research on this problem.at the NACA 
Lewis laboratory, the boundary-layer momentum-integral 
equations are derived and presented herein for a set of 
orthogonal curvilinear coordinates, which may or may not 
be rotating about an arbitrary axis in space and can be laid 
The so gen- 
eralized equations are then transformed by use of an assumed 
velocity distribution and friction law for turbulent boundary 
layer so that an approximate solution can be obtained for 
the boundary-layer momentum thickness and the direction 
of boundary-layer flow. Finally, a numerical solution is 
carried out for the Gruschwitz data in order to make a 
comparison between the estimated and actual measured 


‘values. 


The equations as given in their generalized form are 
readily adaptable to cases of steady, three-dimensional, 
incompressible boundary-layer flow involving centrifugal 
and Coriolis forces. The approximate solution, however, 
has been carried out only for the turbulent boundary layer, 
because in most of the aerodynamic configurations, where 
these equations apply, transition from laminar to turbulent 
flow occurs comparatively early in the flow process. A 
laminar form of the approximate solution can be obtained 
by simple substitution of a suitable velocity profile and 
friction law. , 

It should be noted that whereas the differential equations 
describe the flow phenomena with only the accepted simpli- 
fications, the approximate solution depends to some extent 
on the assumed boundary-layer velocity profiles and the 


relation for friction. Both of these assumptions were made 


1 Supersedes NACA TN 2310, “Generalization of Boundary-Layer Momentum-Integral Equations to Three-Dimensional Flows Including Those of Rotating System” by Artur 


Mager, 1951. 
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on the basis of the data of Gruschwitz (reference 9) only, 
because the data of reference 12 were not adaptable to 
extensive computations for the purpose of this analysis. 
The measurements of Gruschwitz, on the other hand, have 
certain shortcomings as they were obtained in a nonrotating 


channel formed by two circular-arc shaped walls. 


Thus the 


generality of the velocity profiles measured by Gruschwitz 


is in question. 


A revision of the approximate solution can 


therefore be expected when more data become available. In 
addition, any speculation on the occurrence of boundary- 
layer separation (which by definition is a special form of a 
velocity profile) would be absolutely meaningless; no further 
mention will therefore be made of this phase of the problem. 


SYMBOLS 


The following symbols are used in this report (the 
dimensions are given in right-hand column): 


A 


Sl 


Ez) 
Ez) 


Ex) 


rr, rys + 


F 


fx ify „Jz 
G 


constant occurring in second ap- 
proximation for 6 

resultant acceleration vector in 
fixed (inertial) system 

resultant acceleration vector in 
Cartesian coordinate system 

constant occurring in second ap- 
proximation for O 

position vector of particle 


oe 


curvature of z-axis (fig. 1), 5 


constant > ô 
0.01255 Lz 
TELJ) 

z 
(| r 
e\4 Jri 
(E sre y iz) 
e\# 2, U 


rate-of-strain components 

resultant-force vector acting on 
particle 

components of body forces per 
unit mass 

function describing boundary- 
layer velocity profile, also taken 


L 
yy 
as (¥) 


funetion describing boundary- 
layer velocity profile, also taken 


: y A 
as (1 2) 


(=) 
(l) 


d?) 


(0) 
0) 
(mit?) 
(it~?) 


(0) 


(0) 








ICL 


t 
U,V, W 


U, Va Wo 


Sn) S4 f- 
8 


Gr, ba, Ozz, baz 


quantities describing relations 
among various characteristic 
loss thicknesses in boundary 
layer 

transformation coefficients 

length 

parameter determining nature of 
boundary-layer equations 

mass 

static pressure 

components of stress per unit 
area in Cartesian coordinate 
system 

resultant velocity vector 

perpendicular distance of particle 
from axis of rotation 





Reynolds number based on 6,, 


bU 
v 
radius of circle 


total path length f° dz 


arc length 

time 

values of u, 2, and w outside 
boundary lay er 

velocities in Cartesian coordinate 
system 

time averaged velocities in curvi- 
linear coordinate system 

Cartesian coordinate system 

orthogonal curvilinear coordinate 
system 

function used in transformation 

boundary-layer deflection angle 
measured from direction of re- 
sultant skin-friction stress to 
direction of flow outside bound- 
ary layer 

angle between X-axis and tangent 
to x-axis 

boundary-layer thickness 

displacement thicknesses in three- 
dimensional boundary layer 


measure of boundary-layer deflec- _ 


tion, tan a 

slope of characteristic line 

generalized boundary-layer 
momentum-loss thickness, 
ORe 

momentum-loss thicknesses in 
three-dimensional boundary 
layer 


(0). 


(0) 


(0) 


(mit?) 
(mi3) 


(=) 
(D 


(0) 
(8) 
(2 
() 


m’) 
(it~) 
(it~) 
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@) | 





A variable of function 4, 
(2-7 f a iz) 
3 Ti 
y kinematic viscosity (Pr) 
én E components of vorticity vector (i71) 
p density (ml?) 
Try Ta apparent stresses existing in tur- (ml~é-*) 
bulent flow 
To shear stress at wall (mit?) 
Y arbitrary function satisfying 
equation (39) and boundary 
conditions 
GH+9 4 -6E+9) 7 
(eee) 
w angular velocity E) 
Wx, Wy, Wz eomponents of vector w in Car- (é") 
tesian coordinate system 
Oz, Wy, We components of vector w in curvi- (71) 
linear coordinate system 
Subscripts: 
i initial value 
£ z-direction 
Z z-direction 
IL, II order of approximations 


For Gruschwitz data-point designations and streamline 
designations, see figure 2. 
X 





X=-[7 cos £ da + z sin £ 


Y-y, Z=Z, + z cos @ 
Zi =constant-J* sin 8dz 


d 
Ge 88) 


Figure 1,—Transformation from Cartesian coordinates X, Y, Z to orthogonal curvilinear 
coordinates z, y, Z. 


DERIVATION OF BOUNDARY-LAYER MOMENTUM-INTEGRAL 
EQUATIONS 


The equations for steady flow of a fluid having constant 
density are derived in a Cartesian coordinate system X,Y,Z 


272483—54—-14 











along the streamline. 
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(a) Channel and measuriag plate seen from below. 
(b} Measuring plate seen from below, showing point and streamline designations. Shaded 
sections indicate regions of potential-fiow breakdown. 


FIGURE 2,—Experimental setup of Gruschwits (from figs. 1 and 5 of reference 9). 


rotating with uniform angular velocity about an arbitrary 
axis in space. These equations are then transformed to an 
orthogonal curvilinear coordinate system x,y,z such that 
the x-axis can be placed along any convenient path in the 
XZ-plane, which is considered as a plane of a wall. These 
equations are then simplified in a manner consistent with 
the boundary-layer theory. If the path is chosen so as to 
match a streamline of the potential flow, only one velocity 
will exist outside the boundary layer, that along the stream- 
line. Furthermore, the changes in boundary-layer quantities 
in a direction other than that along the streamline are ex- 
pected to be relatively small in comparison with the changes 
Additional simplifications may thus 
be possible. Finally, integration through the boundary 
layer gives the generalized form of momentum-integral 
equations for three-dimensional flows that may or may not 
involve rotation of the system. 

Equations for steady flow of fluid with constant density 
in rotating Cartesian ceordinate system.—The Navier- 
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Stokes equations of flow for a fixed Cartesian coordinate. 
system X, F,Z (reference 13, p. 576) are 

















p oem pfx t Pax 4 Pix | Open (1a) 
p etm pfy-+ bar 4 bre 4 OPey (1b) 
opeta elt ae a 

and the equation of continuity is 
Boe oh ON oo (2) 


It is now assumed that this Cartesian coordinate system is 
rotating with a uniform angular velocity œ and that the 
observations of the motion of the fluid particles are still 
made from a position rigidly attached to the same system. 
The velocity ge and acceleration @, are as seen by the ob- 
server, that is, they are relative to the X,Y,Z system. 
Because of the. rotation, however, the X,Y,Z system is not 
an inertial system (reference 14, p. 53) and thus the second 
law of motion holds only with respect to acceleration a 
relative to some other system that is nonrotating, 


ma=F 
In terms of @, then (reference 14, p. 104), 
miot maX (wXb)+2moXqo=F 


Here mwX(wXb6) represents the centrifugal force and 
2mwX qe is the Coriolis force. 


Thus for a Cartesian coordinate system rotating with a 


f f : DU, DY, 
uniform angular velocity w, the expressions for Te De’ 





DW, ; . 
and py Must be modified by proper components of the 
Coriolis and centrifugal accelerations. For steady flow, the 
component accelerations as referred to a rotating Cartesian 
coordinate system are therefore 





E. E pas OU; nee oR 
T Us 5X T i; oF T FE? 57 +2(wy Fo— wz Vi — aR 5X 
(3a) 

Di OW ON EON ies Or 
Tye Ose Tt Fo oY +h, EWA +2(wzU— wx H wh Sy 
- (3b) 

DW, OV, Ol a Oe oR 





Di => U, srt y o OY +Ws tox Vo— wry U)— oR 
(3c) 


The equation of continuity, which does not involve any 
yccelerations, remains the same. 
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Transformation to orthogonal curvilinear coordinate sys- 
tem.—-Transformations similar to those of Gruschwitz 
(reference 9) are used as indicated in figure 1 with the pre- 
caution that the system remain right-handed. 


x= |" cos Bdx-+z sin 8 
0 


Y=y (4) 
Z=Z;+z cos B 
where 
Z,=constant— if. sin Bdx 
and 


B=B(z) 


Use of these transformations permits an arbitrary curva- 
ture of the z-axis in only one plane, the XZ-plane. Thus 
the solution is somewhat restricted. In two-dimensional 
boundary-layer investigations, however, it is found that the 
boundary-layer equations are unaffected if the radius of 
curvature in the XY-plane is large as compared with the 
boundary-layer thickness (reference 15, p. 120). In three- 
dimensional boundary layer the same limitation will 
probably apply providing, of course, the values of wrs, wy, 
and w, are properly adjusted to take care of this additional 
curvature. Setting 


d x 
=F (curvature of x-axis) 





gives 
oF +ez) cos 8 = dX sin B 
DZ (1+2) sin 6 570 DZ cosg 


The elements of length at (z, y, z) in the direction of the 
increasing coordinates are (reference 15, p. 101): 


hidz, hady, hadz 
Thus, 


(ds)?= (h) (d£) + (hP (dy) + (hs)?(dz)*= (dX)? + (dY + dZ} 


But because 


E eee A 


and so forth, 


(ds}=(1+ez} (dz} + (dy) + (d2} 
and : 


Ay=(1-+e2) hg=1 Ag=1 pa Š (5) 
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The expressions for the linear accelerations ean be written 
directly, as given in reference 13 (p. 158). (It should be 
noted that the A values herein are reciprocals of those in 
reference 13.) The components of a gradient now are 


Ld 
hz Oz 


ro: 
hy oy 


1 3 
h, Ox 





whereas the components of 9 X@ remain 
Vo—We, We—Ue, Uoy— Vor 


Thus the accelerations in the rotating z,y,z system are 
written as 








het aU 4 V (ydh ae 
Dt h dx | Ay Oy | hz hika Dy 
W frr dh i -} ap òE 
pr (U Seo o > aR SS 
And the expressions for n and Pr follow from symmetry. 


The equation for the divergence now has the form 


div B= EEE En, K (hehal) ++ vy ls We aha ) eo 


whereas the components of the curl ge are 


l ò ò 
rlen- an] 


1 d, d > 


1 Ò., ò 
t=rnl se uP —e dal? | 


In order to obtain the viscous terms the preceding expres- 
sions are used in the expansion of 


y [grad (div J.) —curl (eurl Fa) 


If equations (5) are substituted into these general expres- 
sions and the differentiations are carried out, the equations 
for flow in an orthogonal curvilinear coordinate system rotat- 
ing with an angular velocity w are obtained. The body forces 
are neglected here. 














U OU i z 1 pÒ 

I+ pee ths Sots a+ ise" oa 
g 1 1dP oU 

Uw H ON I e p p sr? i Or? 

z DU de DU , vU e òU Uc? 
(1 Fez} òr dz oy? dz? ' I+ez dòz (1+ezř 
OW de 2e on” (6a) 
(i+ez) dz (1+ez)? dz 
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ÒV Ly ye 1 
uss Sets Se Tee 


giae REHU aW) 
ees V z OV ae 
To oe (Fez dr? (i+ez)? ox dz 





























als OV, e dl : 
wot i] (6b) 
if OW dW 0W € 55 4p OR 
l+ze ox oe ape Oz eo ~ana T 
k _10P [ 1 wL 
aWV —wl)=—F 52+" | fea Oz? 
U de z Wade SWW | 
(itezBdzr (itez) ox dr dy? dz” ` 
We ¢ aW 2 a es 
(ltez?!i+tez d2 “(1+ez)? òr 
1 òU, dV dW, We , 
1+ez Dr On +- oz tige” (6d) 


In the general orthogonal coordinates, the expressions for 
the rate-of-strain components are 


1U, V dh, W òh 
hy or hh hiha Oy h ahı dz 
[1d9V, W dh, U dh 
w=? (7; TASATU OF 
og (LOW, U dh, V_ dhs 
a 2 hs Oz Aghy or hohy dy. 
oa 
=% dy Ep Òz he 
eat 2 W pa d ia 
h hidr \h 
ha d 6 
Ch, OL tes h 


The viscous terms in equations (6a), (6b), and (6c) may be 
expressed using the rate-of-strain components as 























dU z. dUde a dU 
á [ary dx? (Fez or dr taat 
EF L- (14 ot aa sex Ky a 
Sheamtatete) ë A 
Ee ao at a batt Te 87 | 
TE ine TE ex | (7b) g 
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| 1 W U de_ z 

$ (l+ez)? òr? (l1tez¥ dz (1+ez2) 
OW OW We e OW Žž 2e a] 
oy? dz? (1+e2}' 1+ez òz (+ez) dz 


1 Dez, , Òeyz , Dez 
Fe epee 


oWde 
or dz 

















ji z (enen) | (7e) 


Equations (6) are directly applicable to the laminar flow. 
For turbulent flow, because of the velocity fluctuations it is 
necessary to modify the stresses by addition of the so-called 
Reynolds’ stresses. Thus, making use of the parallel form 
in equations (7); the Navier-Stokes equations of motion for 
turbulent flow may be written in terms of the apparent 
stresses as 


eau, „ÒU c 
1+ez Erani Sori Oz ti 








l. apo D 
TET” Roz t2% —w,} == a 


I 1 Ocr OT yr OT ex 2¢rx 
p E T n dz He] Ga) 





: REE Rol 
~ 1-+e2 j 


vE+v ws or U—, IV) 

Taart 5 te Orn 

p OY ips Ox 

U OW oo 
I+ez dz ae 








ELELO CT ye | (8b) 


l-+ez 





: c P OR 
Tag ae ae 





1oP Dra y 
2V — ay ua a t an 


OTy , OG, eo 
dy 32 ~T+ez 





(e—a: | (Sc) 


Simplification for flow within boundary layer.—Equa- 
tions (6) and (8) are equivalent to the complete Navier- 
Stokes equations. Within the boundary layer, however, 
certain terms whose contribution is relatively unimportant 
can be neglected. If the y-axis is taken as normal to the 
wall, the boundary-layer flow then takes place over the zz- 
plane (or the XYZ-plane). All terms are now made dimen- 
sionless by referring the lengths to some body length, the 
velocities to their free-stream values, and so forth, as ex- 
plained in reference 16 (p. 45), and all quantities of the order 
of magnitude of 6 or smaller are neglected. Furthermore, 
because the boundary-layer flow along a definite path 
z=0 is of interest, additional simplifications are possible. 
Setting z=0 restricts the equations, because the general 
boundary conditions (not on the z-axis) cannot be satisfied. 
It will subsequently be seen, however, that these general 
boundary conditions are unnecessary in the solution of the 
final equations. These simplifications yield the Navier- 
Stokes equations for flow within the boundary layer in a 





rotating orthogonal curvilinear coordinate system evaluated 
at z=0, 





Ou, OU, Ou POL’, pòk du 
Von Mage ost ee US eee aed Zayw-rv dy? Sat) 
(9a) 
oR ioP 
ai SE = antes 
wh ay T eN wrw) = n Iy (9b) 
ow fà] ow 1 oP Dw 
Usp? p tU = —we= — nO or wre A+ Daur ($8 St) 


(9c) 


for the laminar case. For the turbulent boundary layer, a 
corresponding set of enone is obtained with the substitu- 


2 
Homer O2# fone & =i) and 4 L Dte for 3 sn) 
p "Voy? 


p oy oy? 

Equation (9b) shows, as ee out in reference 10, that 
because all the terms on the left-hand side of the equations 
are of the order of magnitude of one, within the boundary 
layer, P can vary at most by an amount of the order of ô. 
It is reasonable then to neglect this variation and consider 
P solely a function of the flow outside the boundary layer. 
Thus, if æ is chosen to coincide with a streamline of the flow 
outside the boundary layer, V= W=0, and by integration 
of equation (9a) with the effect of viscosity neglected the 
following relation is obtained: 








P=constant—+ pU+5 pak? © (10) 


which is a form of the equation of Bernoulli. 

Furthermore, because outside the boundary layer the flow 
with respect to some nonrotating set of coordinates is irrota- 
tional with reference to the rotating coordinates the com- 
ponents of the vorticity vector become 


E= — 2a, 
n=— 2o (11) 


This assumption of irrotationality is not always true and in 
some applications, such as the later stages of an axial com-. 
pressor, it cannot be used. As long as vorticity is distributed 
according to some definite pattern, however, a relation be- 
tween, the components of vorticity and the components of 
rotational velocity may be found and substituted for equa- 
tions (11). 

Substituting again in the expression for the components 
of vorticity gives 


Peas ee oe 
2a É (ite2)0—2 w] 


And for z=0, W=0, which is along the streamline, the 
expression for curvature becomes 


=- P pe 2) 
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The equation of continuity remains 
(13) 


Generalized boundary-layer momentum-integral equa- 
tions.—In order to obtain the boundary-layer momentum- 
integral equations, equations (9a) and (9c) are integrated 
with respect to y through the boundary layer to some con- 
stant height d such that 


d>6 
ò (C 5 d au d : > pP 
Zf dyt 2o | wdy e f wdy—U Zf E 


d d a 
U 2 f ee f uwdy—4 af uw dy— 
Oz Oz Jo U 


0] 
a tes 


T Pd 
2d T =j x uwdy+2ay, E wdy=Ud = ž (4a) 
Oz JO 0) P 


U 
=f l i f : f, : 0 
Ò p t d j Yy 


d 
These equations apply equally well for the laminar or tur- 
bulent boundary layer, with the value of r, representing the 
shear stress at the wall accordingly adjusted. By suitable 


use of equations (12) and (13), these equations may be 
transformed to 





and 


(14b) 








aN (U— —uudy+Se (U—u)dy+ an (U —u)wdy— 





Jy Ta — 74% P 
5 S4 = {" (U—u)wd — (U—wwdy (15a) 
and 
ò 1 òU 1 òU 
Sk uwdy+ 5z A wdy—T Sz d ‘wy te aah wdy 


y U 





g r 
= dy+2o f" udy+72 T an wily | vay 
(15b) 


The following definitions are now introduced: The mo- 
mentum thickness in the z-direction of the flow in the z- 
direction, 


L pè 
=r f, (U—u)udy (16) 
The displacement thickness in the z-direction, 


E os lial 7 = 
na | U—way an) 


The momentum thickness in the z-direction of the flow in 
the z-direction, 








1 fe, " 
=n, wdy (18) 
The displacement thickness in the z-direction, 


Pee OE i d 
=$ f wdy (19) 


The momentum thickness in the z-direction of the flow in 
the x-direction, 
1 pi ; 
=T3 |, (U—u) wdy (20) 
The momentum thickness in the z-direction of the flow in 
the z-direction, 


i pare 3 
u= |, wuay (21) 


All these thicknesses, as in two-dimensional boundary- 
layer theory, have a dimension of length. Furthermore, 


5 1 g J } d > d 1 d 
=z, | w TA (U —u)w v= |. uwdy =b 
(22) 


With the use of definitions (16) to (21) and equation (22), 
equations (15a) and (15c) reduce, for z=0, to 








06,, 1 dU Ò brz oz 
tra Cet tae —47 0a = aur? (23a) 
and 
00, , OCE — brz) aba) L of ip Se = 
. aas (23b) 


Reduction of equations to forms obtained by other 
investigators. —If only two-dimensional flow exists, that 
is, if e=0, w=0, and w=0, then equation (23b) vanishes and 
equation (23a) becomes an ordinary Kármán momentum- 
integral equation 


ae 1 òU A 
ty Dr Sz Cet [72 





If w=0, that is, if the system is nonrotating, equations 
(23a) and (23b) become identical with the equations of 
Gruschwitz (reference 9). 

Setting c=0 in equations (9a) and (9c) makes these equa- 
tions identical with the equations of Burgers (reference 10), 
who carried out his derivation for a Cartesian coordinate 
system, 

Finally, if the system of coordinates is chosen so as to 


maintain the right-hand rule and c is set equal to 5, thus 


establishing the z-axis as a circle, then a =1 and because 


of axial symmetry all derivatives with respect to x vanish. 
The coordinates are now assumed to be in a fluid that is 
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motionless åt great distance from the surface of the rotating 


immersed disk. Thus, 
U=W=0 wo=0 
Integration of equations (9a) and (9c) gives, after some 
manipulation, 
a ( “umn dy)=—"22 Dart 
dp o p 
and $ 


d ad 


which are identical with equations of Kármán for the 
rotating disk (reference 17). 


APPROXIMATE SOLUTION OF MOMENTUM-INTEGRAL 
EQUATIONS FOR TURBULENT BOUNDARY LAYER 


Transformation and reduction of dependent variables. — 
In order to obtain a solution of the momentum-integral 
equations, additional relations are needed describing the 
velocity profiles existing in the boundary layer and the 
friction at the wall. 

‘With the use of a suggestion by Prandtl (reference 11), 
the expressions for u and w that will be used are 


n= UG (2) 


(24) 
m= lG (4)s(¥)=eug . 
with boundary conditions on @ and g 
for y=, C=}, g=0 
for y=0, Cc=0, g=l 
and with e defined as i 
e=tan a (25) 


where @ is the angle between the direction of the resultant 
skin-friction stress and the direction of the flow outside the 
boundary layer. Because of this definition of e, g=1 at 


y=0 because ' 
ow 
. OY p 
podu” 
oy 
or 


T 0,27 ETo.z (26) 

Mathematically, such use of e implies a linear variation of 
w with e and makes possible the dissociation of the w velocity 
profile from its scale and direction. Because the flow must 
change direction in the boundary layer from that at the wall 
to that in the free stream, there is no reason to assume that 
such a dissociation is actually possible. In other words, 
there is no reason to believe that g should be a function of 
(y/5) aloneand notofeaswell. Inaccordance with reference 11, 
however, this approximation is certainly admissible for 


small values of e and gives results of qualitative accuracy |. 
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for moderately large «. In addition, in order to check this 


: lw i 
assumption, the value of =- for several experimental 


€ L 
velocity profiles and values of e ranging from 0.216 to 
0.670, as obtained from reference 9, is plotted against 4/8 in 
figure 3. The results of this plot indicate indeed that Gg is 
independent of e, l 
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FIGURE 3.—Plot of! a against 4 for various data from reference 9. = 
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In parallel to the two-dimensional boundary-layer theory, 
the following definitions are made: a ae E 


a 
[a-oa 
JO 








) =H 
[a-@@dy 
Q 
(1—6)Ggdy 
Sa =J 
[a-@@dy 
f (27) 
i Ggdy 
0 =K 
[a-@edy 
d 
[i @edy 
Qe Se Se 


f a-oad 


The relations among the various thicknesses may then be 
written 


=H e, 
brz= E SO, 
(28) 
6.—eK6, 
j e 
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The quantities H, J, K, and L are functions of G and g. 
Because Gand g are representative of the velocity profiles in 
the boundary layer, the changes in these velocity profiles 
must be reflected in turn in the values of H, J, K, and L. 
In other words, the external forces acting on the boundary 
layer and influencing the changes in the shape of the velocity 
profiles also cause a variation in H, J, K, and L. Unfortu- 
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— — = Darto of reference l2, position l-@, converted 
to 2, Y, Z system 
.— Data of reference 12, position Ic, converted 
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(a) Experimental velocity profile from reference 9 (point 15, IID. Data obtained in curved 
duet; = 40 millimeters. 
(b} Experimental velocity profile of reference 12 converted to z, y, z coordinate system. Data 
obtained in boundary layer of yawed wing. 


FIGURE 4,—Comparison of assumed G and Gg with experimental velocity profiles. Assumed 
profiles: Gae(y/8)¥?; Gg=(y/8) Y7 (1-y/8)2. 
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o Data of reference 12, position 1-a, 
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Ficguge §5.—Comparison of assumed correction function g(y/8) with data of reference 12 con- 
verted to 7, y, 2 coordinate system. Assumption: g{y/é) = (i-y/8)?. 








nately, 
large changes in the shape of the velocity profiles and the 
quantities H, J, K, and ZL. This behavior of the velocity 
profiles is verified in figure 3. The data of reference 12 do 


indicate large changes in the shape of the velocity profile; 


however, the data are not presented with sufficient detail to 
permit an accurate evaluation of H, J, K, and L. Thus, 
until more extensive experimental data become available, 
the quantities H, J, K, and L are assumed to be constants 


that can be evaluated either by assuming a suitable form for 


G and g or by computing directly from Gruschwitz data. 
In accordance with reference 9, good assumptions for G 


and g are: 
I 
(3Y 
a=(3) 


fN 
g=(1 £) 


An indication of the degree of fit afforded by these expres- 
sions can be obtained from figure 4 (a), where a calculated 
profile with n=7 is compared with one of the profiles of 
Gruschwitz. Other profiles of Gruschwitz data give similar 
results. It should be noted that this good agreement should . 
not be interpreted as meaning that assumptions (29) will 
always give a good representation of the velocity profiles in 
the three-dimensional turbulent boundary layer. Figure 4(b) 
shows a comparison similar to that of figure 4(a) with pro- 
files converted to the x,y,z system using data from reference 
12. Equations (29) do not afford a good fit in figure 4(b), 


although the equations do represent the general behavior _ 


of the velocities. This comparison is further illustrated in 
figure 5, where the value of g (y/5) as obtained by converting 
the profiles of reference 12 to the z,y,z system at indicated 
points is compared with (1—y/8)”. 


With the use of relations (29), H, J, A, and L are com- 








puted as 
H= n 
~ h 
J= n11n+7) 
(Qn+ 1)(8n-+1)(8n+2) 
(30) 
2n?(2+n) 
AS (2n+1)(8n+1) 
6n* 
~ Bn +2)2n+1)(Gn+2) 
which for n=7 give 
H=1.2857 
J=0.5423 
(30a) 
K=2.6727 


£=1.1285 


the available data of reference 9 do not involve 
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Averaging the values along line III of Gruschwitz data 
(fig. 2) results in 


H=1.37 


J=0.550 
(30b) 
K=2.43 


L=0.968 


This relatively good agreement between the two sets of 
values is also indicative of the over-all fit of the assumed 
expressions for G and g to the data of reference 9. 

The additional relation that is needed for the solution of 
the momentum-integral equation is the expression for surface 
friction. In reference 9, Gruschwitz demonstrates that 
Kármán’s friction law 


Foz 
“Hyg 0.01255 (Fz) (31) 


appears to be valid in the three-dimensional boundary layer 
as well. Substituting relations (26), (28), and (31) in 
equations (23) yields i 





06; 6, OU 06, w Toz 
Sa TU SF 57 2+ +de g +0, OS 4 g Jeb ric (82a) 
and 
de 06, L 
21 6,€ set be’ 5 4 52 TE- p oe F oa) <a Jes 
6,0U eoU 2w, oz 
(Lé—1—B) 75+ UK —Jye 5 4 Le’) T b= —e Ter 
(32b) 


Because of the form of the relation for friction, an advan- 
tageous transformation of variable is 


O=0, Ret (33) 


in order to eliminate the Reynolds number from the equa- 
tions. 

With the use of equations (33) and (12), two nonlinear 
partial differential equations for O and ¢ applying along z=0 
are obtained from equations (32), 





400, 4 5H+9\ 1 dU 
Zor 5 Jes Seto Si AG los 
Ti Sye 8—0.01255>=0 (34a) 
5 U 
and 
oe 4 


5 K- Neste Le 24K Te S210 È sot 


a Ke AE oF TEE Le -1—H) io 





2-H) 72040. 012556 |=0 (34b) 











As shown in the appendix, these equations can be either 
hyperbolic, parabolic, or elliptic, depending on the shape of 
the velocity profiles existing in the boundary layer. For 
u=U(y/6y"" and g=(1—y/6)?, the equations are always 
elliptic. 

Simplification of equations and approximate solution.— 
The relative importance of the various derivatives in equa- 
tions (34) is now determined. First, O and «e are assumed to 
be quantities that are smaller than one, which can be accom- 
plished simply by referring all lengths to a total path length 
S and a to 45°. As a result of this assumption, all deriva- 
tives in © and «e become of the order of magnitude of O or e. 
Equation (34a) is then divided by 4/5, so that the coefficient 
of 00/dz is 1 and the coefficients of 00/Oz and dedz arc Je 


and 5 Jo, that is, of the order of magnitude of e and ©, 


respectively. In a similar way equation (34b) is divided by 
(K—JO so thatthe coefficients of 00/dz and dedz become 


S e/O and 1 (order of magnitude of one), respectively, and 


2 
those of 06/Oz and d¢/dz2 become $ L 5 and 2L e (order of 


magnitude of e), respectively. Then, if e is small as com- 
pared with tan 45° and © is small when compared with S, all 
terms of the order of magnitude of ê, Oe, and O? may be 
neglected, which gives 





als) 5H+9\ 1 òU] 
and 
Oe 360850 1 ÐU , 0.01255 1 
o¢'|5002'5U dz a J Od 
is 2 wy ; 
=J oop [Othe Ho (35b) 


These two expressions show that the primary changes in 
© and e in equations (34a) and (84b) occur only in the 
x-direction and thus the description of the phenomena only 
z=0 is justified. 
A solution of equations (35) can now be obtained by suc- 
cessive approximations because U, wy, and c are assumed to 
be known functions of x. First, equation (35a) is solved, 


un) , 0.01569 fe (=) 
T) tay [Uo 7de 0 
ae tae 


The values of ©;(z) are then used in an appr oximate solution 
of equation (35b) 


O;=0; 


4 


9 
o [UN (OY 
ey= é; Ə: 


T) Bet i 


(J—KE) 6U PE 1 (2) 


z £8 
fomen e end en 


where 





0.01255 
E (e= eE- SE a - ` (38) 
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With e(z) known, the approximation for O could in turn 
be improved by again solving equation (35a). If œ grows 
large along the path, however, it is more advantageous to 
consider the following equation: 


d, , 00,7 /5H+9\1 dU Ta 
a (a sete 





-tigen 8=0.01569 


(39) 


U 


It is thus hoped that neglection of de/dz will notaffect theac- 
curacy of the solution to a very marked extent. The solution 
of equation (39), which may be obtained by the method of 
Lagrange, is 


On=| vf = edz) | no : 4 





(== 2) 
0. aal f US * / Rila) Ti 


z (40) 
gaga PO. 
where 
z = Cedr i 
Toar (41) 
18 w 
qd | eas 
E (2) =e* f (42) 
and 


W(2-J f eaz) 


is an arbitrary function satisfying the boundary condition; 
when =r; and z=0, then On=9;. Setting 


(2-3 [edz )=a 


atz=x;,givesh\=0. In addition, for z=0, dis of the order of 
magnitude of e. 


Expansion of y in Maclaurin’s series about X(x;) yields 
rw 
Y= (0) 0z y” 0) 


Because there is only one boundary condition, it is possible 
to determine only one of the constants in this expansion; 
consequently, y cannot be uniquely established. The fact 
that à is of the order of magnitude of e, however, suggests 
that the assumptions made for 4” (0), ¥’’’(0), and so forth, 
are successively less important. Thus these derivatives may 
arbitrarily be expressed by a single constant, 


oe a E), a 


(43) 


H+9 
¥O)= =o (eaa 


where A from purely dimensional considerations must have 
the dimensions of /-!, From expressions (41) and (42), it is 


suspected that 
A=B(5-3 a) (44) 
i 





where B must be obtained from the experimental measure- 
ments. 

It should be noted that because ¥(A) cannot be uniquely 
determined other functions of \ satisfying the single boundary 
condition could be used as well. The function e4* is chosen 
only because it is convenient to use and parallels the expres- 
sions (41) and (42). This arbitrariness of the functional form 
of y and the value of A is due to the consequences of assuming 
z=0, and thus it is probably not advisable to carry any 
further approximation for eg and so forth. 


In solving equations (36) to (44), either set of values for _ 


H, J, K, and L may be used. Because the averaged values 
(30b) are probably more representative, having been obtained 
by evaluating experimental data at a number of different 
positions, it is advantageous to use these values in com- 
putations. 


COMPARISON WITH EXPERIMENT 


In order to check the validity of the approximate solution, 


the boundary layer along four streamlines of reference 9 was __ 


computed and compared with the measured values. The 
designation of the streamlines and data points is illustrated, 
in figure 2. Because the data were taken along curves I to 
V of figure 2, the computation along a streamline requires 
first an interpolation among the various data points. As a 
result of this interpolation, the computations could not be 
carried through the full length of each streamline. Values 
(30b) were used for quantities H, J, K,and L. The constant 
B was obtained by fitting along streamline B the solution 
for Oy, so that at z=S, Ou~O measured. In this manner, 
the value of B was found to be 38.5. This value was then 
used in computations of streamlines A, C, and D. It is 
noted that B=38.5~7(Re™*),, although justification for such 
a dependence cannot be made. In all integrations Simpson’s 
rule was used. 


The results of the computations are plotted in a nondi- 
mensional form and compared with the interpolated measured 
values in figures 6 and 7. A study of these figures reveals 
a fair quantitative agreement between the measured and 
estimated values of O and a. As the values of e-tan 45° 
(fig. 7) the first approximation for © in figure 6 becomes 
progressively worse, which is remedied by the second approx- 
imation. The poorest agreement is obtained along stream- 
lines A and D, which because they are closest to the walls 
might be affected by the flow in the corners of the duct. 
Streamline D especially may be affected inasmuch as 
Gruschwitz mentions the existence of separation on the 
convex wall. 


The fair quantitative agreement with the measured values 
is not to be interpreted as a conclusive check of the validity 
of the procedure and the assumed values in all cases of three- 
dimensional boundary-layer flow. The suggested procedure 
simply represents the best that can be done in view of the _ 
meagerness of the available data. Because the Gruschwitz 


data do not involve the effects of uniform angular velocity 
1 òU 1 dU 


and because the variations in = and = — are small, 


U dr U dz 
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Distance al ong streamline, xS 


(a) Streamline A; S=22.8 inches; 6;=0,282 inch. 
(b) Streamline B; S=34.75 inches; 6;=06,289 inch, 
(c) Streamline C; S=37,8 inches; 6;=0.286 inch. 
(d) Streamline D; S=21.25 inches; 6;=0.244 inch. 
Fieure 6.—Comparison of calculated and measured generalized momentum thickness. 


Experimental data from reference 9. 
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it could be maintained that this check of the procedure has 
been carried out on a somewhat special case. For that 
reason, it is desirable that additional experiments be carried 
out in setups that eliminate the present shortcomings. A 
larger variation of Reynolds number should also be used. 
With additional experiments, a modification of the values 
of H, J, K, L, and B, together perhaps with some refinements 
of the procedure, will be in order. It might be well to remem- 
ber, at such time, that because of the necessary empiricism 
involved (which results from the very limited knowledge of 
turbulent phenomena), long and tetlious computations would 
rarely be worthwhile. 


CONCLUSIONS 


The following conclusions can be drawn from an analysis 
of the three-dimensional momentum-integral equations and 
a comparison of the numerical results with the Gruschwitz 
data for turbulent. boundary layer: 

1. Within the boundary layer the static pressure can vary 
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(a) Streamline A; S=22.8 inches. 


at most by an amount of the order of magnitude of the 
boundary-layer thickness 6. 

2. It is possible to generalize the velocities in the boundary 
layer by use of two characteristic quantities ô and e where 
e is the tangent of the angle enclosed by the direction of the 
resultant skin-friction stress and the direction of the flow 
outside the boundary layer. ; x 

3. When the generalized boundary-layer momentum-loss 
thiekness © is small as compared with the total path length 
and e is small as compared with tan 45°, the primary changes 
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- (b) Streamline B; S=34.75 inches. 


Figure 7.—Comparison of calculated and measured bound ary-layer deflection at wall. Experimental dats from refercnee 9, 
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in © and e occur along the streamline of the flow outside the 
boundary layer. 

4, The three-dimensional boundary-layer momentum- 
integral equations can be either hyperbolic, parabolic, or 
elliptic, depending on the relative magnitude of the parameter 
MN, which in turn depends on the shape of the velocity 
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(c} Streamline C; S=37.8 inches. 
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profiles existing in the boundary layer. The power-law 
profile when used with the correction function g=(1—y/8)? 
always results in elliptic equations. 

5. The approximate solution of the three-dimensional 
momentum-integral equations shows a fair quantitative 


agreement with the values measured by Gruschwitz. 
6. Additional experimental data are necessary to establish 


More generally applicable values for form parameters HM, 


J, K, and L and B, the constant used in the second approxi- 
mation for 9. 


Lewis FLIGHT PROPULSION LABORATORY 
NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
CLEVELAND, Onto, November 1, 1950 





32 D or a 
/ 
28 
——— Measured / 
— —~~—Calculated 4 T 











X 
~ 
N 
K ‘N 


S 















~ 
N 





Boundary-~layer deflection, æ, deg 
r ~ 
D 




















F 


l / 
A 


OET 


8 2 rod 40 
Dee along See 2fS 


























(d) Streamline D; S=21.25 inches, 


FIGURE 7, Concluded.—Comparison of calculated and measured boundary-layer deflection at wall. Experimental data from reference 9. 





APPENDIX 


ADDITIONAL REMARKS ON THREE-DIMENSIONAL 
BOUNDARY-LAYER MOMENTUM-INTEGRAL EQUATIONS 


In order to obtain the approximate solution of equation (34), 
it was shown by comparing the relative order of mag- 
nitude of the coefficients that some of the terms may be 
neglected. Care must be taken with such simplifications 
inasmuch as various implications of the equations in question 
may be obscured by this procedure. For this reason, aside 
from the approximate solution, the character of equations (84) 
was also investigated in detail. 

By use of the procedure outlined in reference 18 (p. 38), 
along z=0 the system of equations (34) is found to be 
hyperbolic when MND1, elliptic when MN<1, and para- 
bolic when AfN=1, where 


L AA 


a (K ve J) I bz Ozz 


=0 (45) 





Because J, K, and L are functions of G and g, the character 
of equation (34) depends on the shape of the velocity profiles 
in the boundary layer. 

It should be noted that when MN =0, then L=0, which is 
only possible if G=0 or g=0, and in turn u=0 or w=0. If 
the trivial case u=0 is neglected, it is established that when 
w=0, e=0 as well. But for e=0 and w=0, equations (34) 
reduce to a special case 





pee eu set) 0.01255 (46a) 
5o0r  U Oz 
and 
2w, 1+H 0.01255 
r pg o 5 (46b) 


Here equation (46a) is an ordinary two-dimensional boundary- 
layer momentum-integral equation for © and equation (46b) 
is a relation that evidently must exist among U, 0U/oz, 
wy, and O, when e=0 and w=0. 

When the equations are elliptic, no real characteristic 
direction can be found. When only one characteristic direc- 
tion exists, the equations are parabolic and in the hyperbolic 
case two characteristic directions through each point of the 
zz-plane are obtained. For the parabolic case then, 


d 
Tohe 


and for the hyperbolic case, 


dz L+yit—JK—SL 
_ (K—J) 

















and the characteristic lines are asymmetric with respect to 
the z-axis. In order to determine whether elliptic, parabolic, 
or hyperbolic equations apply, the magnitude of ALN is com- 
puted. Substituting from expressions (30), ALN is obtained 
in terms of n: 

6(3n-+1)(8n+2) _ 1 108n+96 








MN = ntin A t 55 555n*+ 57414) 


This equation shows MAN to be a monotonically increasing 
function ofn. Forn=0, 


MN= 


“1 


and 

een E 

leo 55 
These results indicate that a so-called power-law profile when 
used with g=(l-y/é)? always results in equations that 
although elliptic are very near to being parabolic. Using 
values (80b), 

ALN=0.936 


which again indicates an elliptic character of the equations. 
It should be remembered, however, that the assumptions for 
G and g were made on the basis of only one set of data; con- 


sequently there is no assurance that the velocity distributions j 


existing in the boundary layer will always give the same 
values of MN. In fact, it is generally more likely that they 
will not give the same values of MN. 
variation of MN may already be obtained from figure 8, 
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Fiacre 8.—Values of parameter MN for data of reference 9. 
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Some indication of the _ 
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where its value was plotted for each Gruschwitz data point. 
At points 5, 6, and 7, however, small values of w may have 
prevented an accurate determination of Z and as a ĉon- 
sequence AfN=0 there. The value of MN in figure 8 
varies within the limits 0.65< AfZN< 1.2, with the bulk of the 
points indicating that MN ~0.95. 

On the basis of the preceding discussion, there is some 
evidence of the equations being parabolic, elliptic, and hyper- 
bolic in the turbulent boundary layer. It is interesting to 
note that generally (as in supersonic and subsonic flow, for 
instance) these hyperbolic and elliptic regions have their 
counterpart in physical phenomena. Thus some essential 
differences might exist in the process of momentum transfer 
between the hyperbolic and elliptic regions.. These differ- 
ences cannot now be ascertained because first equations sim- 
ilar to (34) with z0 would have to be obtained, and there 
is no mention of any irregularities in the behavior of the flow 
in reference 9. When additional experiments are made, 
however, it would seem advisable to study closely these two 
mathematical regions in order to obtain some indication of 
the physical make-up of their differences. . 
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